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f-H . Abstract 

We show that the local equivalence problem for second-order ordinary differential 
equations under point transformations is completely characterized by differential invari- 
ants of order at most 10 and that this upper bound is sharp. We also show that, modulo 
Cartan duality and point transformations, the Painleve-I equation can be characterized 
J> ' as the simplest second-order ODE belonging to the class of equations requiring 10th 

order jets for their classification. 

1 Introduction 

o; 

■ Two second-order ordinary differential equations 

q = f(x,u,p), Q = F(X,U,P), p = u x ,q = u xx , P = U X ,Q = U X x 
^ ■ are (locally) point equivalent if there exists a local diffeomorphism 

S; (X,U) =ijj(x,u), ^eDiff(]R 2 ), (1.1) 

such that 

F(X, U, P) = Q(p, f(x, u,p),X x ,X u , U x , U U ,X XX ,X XU ,X UU , U xx , U xu , U uu ) (1.2) 
where 

D x = d x + pd u + qd p 

is the usual total derivative, and 

D X U pU u + U x 



P = P(p,X x ,X u ,U x ,U u )- 

D X X pX u + X x (1.3) 

Q — Qi.Pl Ql -X-Xl -X-ui U x , Uui X-XXl X-XUl X-UUl U xx , U X ui U uu ) 
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D X P _ D Z X U ■ D X X - D X U ■ D 2 X X 



D X X {D x Xf 

{pX u + X x ) (p 2 U uu + 2pU ux + U xx ) - {pU u + U x ) (p 2 X uu + 2pX ux + X xx ) 



+ 



(pX u + X x f 

U U X X U X X U 
( P X U + X X )3 q 



describe the transformation law for u x and u xx , respectively. 

Given two smoothly defined ODEs q = f(x,u,p) and Q = F(X,U,P), equations 
(|1.3p together with the obvious contact conditions, constitute a quasi-linear system of 
PDEs in the group jet variables X, U, X X ,X U , U X ,U U - This system is over-determined, 
owing to higher order integrability conditions that, for sufficiently high order, govern 
the outcome of the equivalence problem. 

Taking repeated derivatives of (jl.2p with respect to X, U, P yields the following 
necessary conditions for equivalence: 

d^F(X U P) ^ 

Fijk(x, u, P) = QxiQjjjgpk = Qijkip-, ix,qu, q P , ■ ■ ■ ; x x ,x u , u x , u u , . . .), 

where 

Q l]k { P ^ +k) ;4 +3+k+2 \ e Diff(M 2 ), 

describes the transformation law for the partial derivatives under point transfor- 
mations (jl.ip . and where 



Uijk ■ < i + j + k < n}, q ijk := q x i^ v k ; 
Ly, C/fj) : < i + j < n}, Xij := X x , uj , L 4J 

denote the indicated jets. 



V> (n) = {(XijtUij): 0<i + j< n}, X i:j := X x i uj , U l3 := U x , u r, (1.4) 
4 n) ={(X ij ,U ij ): l<i+j<n}; 



Definition 1.1. We say that a second-order ODE q = f(x,u,p) is classified by nth 
order jets if the algebraic consistency of the system 

P = P(p,^ ] ), 

F ijk (X, U, P) = Q ijk (p, f^ +k \x, u,p),^ +J+k+2) ), 0<i+j + k<n; 

is sufficient for the existence of a point transformation relating q = f(x,u,p) and 
Q = F(X, U, P). We call the smallest such n the IC (invariant classification) order of 
the differential equation. 

Note: In the formulation of the above definition it must be understood that the pseudo- 
group variables (X X ,X U , U X ,U U , . . .) are to be treated as auxiliary independent variables 
rather than functions of x and u. 

The equivalence problem for second-order ordinary differential equations under point 
transformations was first considered by Lie, Tresse and Cartan in [U [T2| 124] . and was 
revisited by many subsequent authors; see |15l Chapter 12] for a modern treatment and 
for a recent application. The question that motivates us here is the following: 

What is the maximal jet order required for the invariant classification of a 
second-order ordinary differential equation up to local point transformations? 
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By way of an example, the well-known Linearization Theorem for second-order 
ordinary differential equations [H [TOj ESI [22] states that q = f(x,u,p) is point equivalent 
to the trivial equation Q = if and only if the fourth-order relative invariants 

Qpppp = 0' ,., 
^ 9 ^ ^ (1.5) 

D x(Qpp) ~ 4D x (q up ) - q p D x (q pp ) + 6q uu - 3q u q pp + Aq p q up = 0, 

are identically zero. From this it follows that the linearizable class has IC order equal 
to 4. On the other hand, the invariant classification of general second-order equations 
will require higher order jets. The complete answer regarding the maximal order is 
given below in Theorem 11.31 

A well-posed equivalence problem requires some notion of regularity. Therefore, 
before proceeding further, we must impose a technical rank assumption. Owing to the 
covariance of the transformation laws, the functions Qijk are invariant with respect to 
point transformations. For a smoothly defined ODE q = f(x,u,p), let 

q^=f^(x,u,p) 

denote the nth order jet of the defining function. The composition of Qijk and 
produces an invariant of the ODE, 

Qijk = Qijk (p,f^ +k) (x,u,p),4 l ^ +k+2) ) , (1.6) 

which we call a lifted invariant to signify its dependence on the pseudo-group variables 
V'o ^ = (x x ,x u , U x , U u , . . .). 

The above remarks lead us to the following definition, after which we will be ready 
to state our main result. 

Definition 1.2. We say that a smoothly defined ODE q = f(x,u,p) is regular if for 
every n > the rank of the lifted invariants {Qijk '■ 0<i + j + k<n} is constant. 

Theorem 1.3. Every regular second-order ODE is classified, relative to point trans- 
formations, by its 10th order jets. This bound is sharp — meaning that there exist 
regular ODEs that are not classified by 9th order jets. Furthermore, every ODE having 
the maximal IC order of 10 is equivalent, modulo point transformations and Cartan 
duality (see Appendix El for the definition), to an equation of the form 

u xx = 6u 2 + g(x) with D xx [g(x)- 1 / 4 ] ^ 0. (1.7) 

To prove Theorem ll.3| we use the method of equivariant moving frames, |19] , which 
is a theoretical reformulation of Cartan's moving frame theory. The proof boils down 
to identifying the branch(es) of the equivalence problem where non-constant absolute 
invariants appear as late and slowly as possible during the course of Cartan's normal- 
ization procedure. Taking advantage of the universal recurrence relations, we first give 
a symbolic proof of Theorem II .31 Using the fact that there is a notion of duality among 
second-order ordinary differential equations (see Appendix [A]) our conclusion is that 
there exist two families of differential equations (dual to each other) that achieve the 

contrast, absolute differential invariants are functions of x,u,p only. To construct absolute invariants 
one eliminates the group variables from the lifted invariants of sufficiently high order. 
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maximal IC order. By a generalization of Cartan's Integration Theorem, [3j, we are 
then able to show that one of the two families of differential equations depends on 
one arbitrary function of the independent variable. Finally, we integrate the structure 
equations for the canonical coframe and derive form (jl.7p . 

Second-order ordinary differential equations equivalent to ()1.7|) by point transfor- 
mations and Cartan duality admit three fundamental absolute invariants Ij,Is,Ig of 
the indicated order, and a tenth order invariant I±o which is functionally dependent on 
Ig. Relative to the normal form ()1.T|) . these invariants can be expressed as 

T P j 9(x) T _ {g'{x)) A _ g(x)g"(x) 

The functional relation between Ig and I\o is the essential classifying relation for equa- 
tions of maximal IC order. 

The class of equations (jl.7p requiring 10th order jets for their classification includes 
the Painleve-I equation as the subclass when g(x) is linear in x. An invariant char- 
acterization of the Paineleve-I and II equations up to fibre preserving transformations 
was given in [9]. The invariant characterization under the larger pseudo-group of point 
transformations was derived in [11]; see [1] for further discussion. As a Corollary to 
Theorem 11.31 we are able to give the following simpler characterization of Painleve-I. 

Theorem 1.4. The equivalence class of the Painleve-I equation can be characterized 
as the subclass of second-order ODEs requiring 10th order jets for their classification 
and satisfying 

Qpppp = 0, ho = 0. 

The condition q vvvv = distinguishes the Painleve-I equation from it's Cartan dual. 
The vanishing of I\q = means that the Painleve-I equation can be characterized as 
the "simplest" second-order differential equation belonging to the class of maximal IC 
order equations. We will derive (\1.8\i at the end of Section 14.21 . Thereafter Theorem 
11.41 follows as a straight-forward Corollary of Theorem 11.31 



2 Preliminaries 

2.1 The geometric setting 

In this section, we introduce the geometric setting for the equivalence problem of second- 
order ordinary differential equations. The key formalism is a certain groupoid and two 
sets of fundamental equations: the universal recurrence relations for the prolonged jet 
coordinates, and the Maurer-Cartan structure equations of the infinite-dimensional Lie 
pseudo-group V = Diff(R 2 ) [321 [TO]. 

Let N 3 = J^R) and M 4 = J 2 (R,M) denote, respectively, the first- and second- 
order jet space of curves u = u(x). We identify a smoothly defined second-order or- 
dinary differential equation q = f(x,u,p) with a smooth section /: iV 3 — > M 4 and 
denote by = J n r(A^ 3 ,M 4 ) the bundle of nth-order jets of second-order ODEs. 

The corresponding prolongation of the differential equation q = f(x, u,p) is denoted by 
f(n). N 3 _^ j{n)_ Next] let g c Di ff( R 4) denote the prolongation of Diff(lR 2 ) to M 4 . 
For < n < oo, let 2?( n ) — y M? denote the bundle of nth order diffeomorphism jets of 
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Diff(M 2 ), and let — > J^> denote the nth-order lifted bundle obtained by pulling 
back £>(™+ 2 ) £?(») -> ]R 2 via the projection J"( n ) -> M 4 -> iV 3 -> M 2 . 

The jets of g and (X, J7) = ^(x, u) give canonical bundle coordinates as follows: 

• iV 3 : x, u,p; 

• M 4 = J° : x,u,p,q; 

• J** 1 ' : x,u,p,q^; 

• : x,u,X,U,^ ] ; 

• : x,n,p,g( n ),X,J7,4 n+2) - 

We note that the nth-order diffeomorphism jet = {X, U, X x , X u , U x , U u , . . .) 

has 

2 x (1 + 2 + 3 + . . . + n + 1) = (n + l)(n + 2) 

components and that the nth-order jet = (q,q x ,Qu,Qp, ■ ■ •) has (™3 3 ) components. 
Hence, 

dim = 3 + f" g 3 J , dimpW = 2 + (n + l)(n + 2), 

dimfW = 3 + J" ^ + (n + 3)(n + 4). 

Next, we define the groupoid structure of the lifted bundle fW. The source map 
(j( n ) ■ £( n ) j{ n ) is the standard projection given by (x,u,p,q^); while the target 
map T^ 1 ' : £( n ) — > J n is the projection given by the prolonged action (X,U,P,Q( n ty, 
where 

P = 

as per fjl .3|) and 

Q<»> = Q^\p,q^\^ +2) ) = [Q l]k {p,q^\4 +1+k+2) ) :0<i+j + k<n}. 

The point transformation pseudo-group Diff(R 2 ) has a prolonged action on 
and two dual prolonged actions on £^ n \ Given a point transformation (X, U) = f(x, u) 
we define the prolonged actions 

^n) . j(n) ^ j(n) ^ ^(n) . g („) ^ f („) ^ ^(n) . f(n) ^ f(n ) 

according to [18J : 

: (s, u,p, g (n) ) ^ (X, [/, P(p, QW (p, ?W <^ n+2) )) ; 
<^ n) : (x, U ,p,gW,^)^(x,n,p,gW,( ¥ ,oV) {n) ); (2.1) 
: (*, u,p, gW, ^)) ^ (X, 17, P(p, 4 X) ), QW(p, gW ^ n+2) ), o v?" 1 ) <">) . 

From the above definitions, it follows immediately that the source projection tr^ is 
(^-invariant and c^R-equivariant: 

ffW o ^ = ff W o ^ = o ; 

and that, dually, the target projection is <^R-invariant and <^£,-equivariant: 

r (") o <p<g> = T W, r< n » o ^ = OT W. (2.2) 
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2.2 The invariant coframe 

A coframe on J7~(°°) is given by the basic horizontal and contact one- forms 
dx, du, dp, 9 ijk = dq ijk - q i+ i,j,k dx - q { >j+1)k du - q^k+i dp, i, j, k > 0. (2.3) 

The standard coframe on is spanned by dx, du together with the group forms 

Tjj = (7.V, ; A',.;., dx-Xij+i du, = dUij-U i+ ij dx-U itj+ i du, i,j > 0, (2.4) 

while a right-invariant coframe for consists of the one-forms 

w x = X x dx + X u (in, = f7 x fix + U u du (2.5) 

and the Maurer-Cartan one-forms 

The latter are defined, implicitly, by taking formal derivatives of the relations 

T = p, * = x/, 

with respect to x, u, and then solving for the pij, Vij. The first few relations that result 
are shown below: 

T x = X x (jl x + U X uu, = X x v x + U x uu, 

T u = X u p x + U u fiu, 9 U = X u u x + U u vu, 
^Ix = X x pxx + 2X X U X pux + U x puu + X xx px + U xx pu, 
^ux = X X X U pxx + (X X U U + X u U x )pux + U X U U puu + X ux px + U ux pu, 
^uu = -^u Vxx + 2X U U U pux + U u puu + X uu px + U uu pu, 

with the coefficients of the higher-order relations given by the multi-variate Faa-di- 
Bruno polynomials [21] . The precise coordinate expression for the Maurer-Cartan 
forms is derived in |18j . but these are not necessary for the symbolic implementation of 
the moving frame method. 

To construct an invariant coframe on we note that the space of differential 

forms ft* = f2*(£(°°)) on the infinite-order lifted bundle decomposes into 

«* = ® k ,i n k ' 1 , 

where k indicates the number of jet forms fj2.3j) and I the number of group forms (|2.4[) . 
Let fij = 0^ Q k, ° denote the subspace of jet forms, and define the projection map 
7Tj : ft* — > n*j onto the jet component. 

Definition 2.1. The lift transformation A: ft*(J°°) — > fi} is defined by 

A = 7t j o(t( °))*. (2.7) 
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Thus, by construction, the lift of a jet form on J^^ 00 ) is an invariant jet form defined on 
In particular, we have 

u> x = X(dx), uj u = X(du). 
We also introduce the invariant one-forms 

UJ P = \(dp) = TTj(dP) 

= ( P x u + x x )- 2 [(x x u u - X u U x )dp+ 

+ {p 2 {U ux X u — X UX U U ) + p(U xx X u — X XX U U + U UX X X — X UX U X )\ + 
+ {U XX X X — X xx U x ) S jdx ^2 
+ [p 2 {U uu X u — X UU U U ) + p(U uu X x — X UU U X + U UX X U — XuxUu)^ + 
+ (U UX X X - X ux U x )jduj, 

^0,0,0 = A (6*0,0,0 ) = ^y" T ( dq ~ qxdx ~ qudu ~ q P dp ^ > 

and more generally, 

$ijk = H@ijk)- 

Next, we introduce the infinitesimal generator 

v = Z\x,u)—+f(x,u)—+f(x,u,p)— + <P{x,u,p,q)—+ <*> k 



dx du dp dq i+ f^ k ^ dq ijk 

d d d 

v = u)— + rj(x, u )-q^ + [Vx+p{Vu ~ £ x ) ~ P 2 ^q^ + fe^ + q ( 7 lu ~ 2 6<0 

+ p(2% u - Cxx) - 3pgC u + P 2 (r]uu ~ ^ xu ) - AuJtt- + Yl ^ 



dq i+ £z>r ^* 

(2.9) 

of the Diff(K 2 ) action on obtained by prolonging (|1.2j) . The coefficients <^ jffc are 

defined recursively by the usual prolongation formula 

4> ijk = DlDlDfa - £V - fq u - fqp) + ^q ijtk+1 + ^ij+l.fc + fq^k+l, (2.10) 

where 

are the total derivative operators on J'( co \ The definition of the lift map (|2.7|) can now 
be extended to the vector field jets by setting 

■M^ieW) := t^ij-i ^■{ r lx i ui) := v iji hj — 0, 

where /^j, fjj are the Maurer-Cartan forms defined above. We also note that the lift 
of the source variables gives the target variables: 

X = X(x), U = X(u), P = X(p), Q lJ k = X{q ij k). (2.11) 
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Proposition 2.2. The universal recurrence relations for the lifted invariants are 
dX = u) x + /j,, 

oIU = uj u + v, (2.12a) 

dP = UJ P + V X + P(l>£/ - (ix) ~ P 2 VU, 

dQijk = Qi+i,jk" x + Qi,j+i,ku u + Qi j)k +iu p + %fe + A(</> ij ' fc ), (2.12b) 
where \(4>^ k ) is the lift of the prolonged vector field coefficient (|2,10p . 

Proposition 2.3. The Maurer-Cartan structure equations of the diffeomorphism pseudo- 
group (jl.ip are, [20] , 



(0,0)<(M)<(i,j) 



(2.13) 



(0,0)<(fc,€)<(i,j) 



- Aw 1 - A u u . 

3 The equivalence problem 

3.1 The direct, infinite-dimensional formulation 

Given a smoothly defined ODE q = f(x,u,p), let £^ — > iV 3 denote the bundle over 
N 3 given by the pullback of £^ by . The canonical bundle coordinates on £^ are 
x, u,p, X,U,ip^ +2) , which means that 

dim£:j n) = 3 + (n + 3)(n + 4). 
The corresponding embedding £y <^-> is given by 

(x, u, p, X, U, 4 n+2) )^(x, U ,p,f<- n \x,u,p), X, U, 4 n+2) ) . 

(n) 

It follows that £j is no longer a groupoid, but merely a fibre bundle. The correspond- 
ing restrictions of the source and target projections to £f, denoted Of : 8j — > N 3 
and tJ 1 : 8 J 1 — > J^ n \ respectively, are given by 

W . L „ „ v rr „/>+ 2 )> 



ay : (x,u,p,X,U,% ')>-> (x,u,p), 
rf: (x,u,p,X,U,4 n+2) ) ^ (X,U,P(p,4 l) ),Q^(pJ^(x,u,p),^ +2) )). 

Since the mapping 
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has constant rank, the rank of r \ is constant if and only if the ODE is regular as per 

Definition 11.21 Also, by (12.2|) . the restricted target map, is invariant with respect 
to point transformations. Hence, if a given ODE is classified by nth order jets, then 
the imgTj C serves as a signature manifold for the ODE; that is, two ODEs are 
locally point-equivalent if and only if their signatures overlap on an open set. 

The contact invariant one-forms in (|2.12bp are essential when working with the 
invariant variational bicomplex, [3 [23], or studying geometric submanifold flows, [13} 
116] , Since the restriction of £^ to £^ annihilates the one-forms fiijk, the equivalence 
problem must be formulated in terms of the invariant -coframe 

u; x ,u; u ,u} p ,fj,ij,u i j. (3.1) 

The structure equations for this coframe consist of (I2.13P as well as 

duj x = —dpL = fj,x Aw^/ij/A uj u , 
du u = -dp = v x A uf + vn A uj u , 

(3 2) 

cLjP = {v v -[i x - 2Pfiu) Aw p | {u ux + P{yuu ~ Vux) - P 2 ^uu) A u u 
+ iyxx + P(yux - fJ,xx) ~ P 2 Hux) A u x . 

The latter are obtained by computing the exterior derivative of (|2.12ap . and taking into 
account the Maurer-Cartan structure equations (j2. 13|) . 

On the universal recurrence relations (|2.12bp must be considered modulo the 

contact one- forms # = {$ijk}- These relations then express the one- forms dQijt as 
invariant linear combinations of the coframe (|3.ip . So in effect, the present setting can 
be considered as an infinite-dimensional overdetermined equivalence problem; see |15} 
p. 297] for a discussion. 



3.2 The universal reduction 

The direct formulation of the equivalence problem outlined in the preceding subsection 
suffers from an essential difficulty stemming from the fact that the lifted invariants 
depend on an unbounded number of pseudo-group variables. Indeed, there is no upper 

(n) 

bound on rankr^ as n — > oo, and so, apriori, it is not even possible to assert that an 
IC bound exists. 

To overcome this difficulty, we introduce a partial moving frame j!7l [25] for the 
action of the point transformation pseudo-group on In effect, this moves the 

equivalence problem from an infinite-dimensional setting to an 8-dimensional principal 
bundle. After this universal reduction, which is valid for all smoothly defined ODEs, 
the invariant classification proceeds using the usual method of reduction of structure, 
albeit with a certain amount of branching. 

To illustrate the normalization procedure, we first consider the order zero normal- 
ization in some details, and then pass to the description of the full normalization. 
Considering the recurrence relations (|2.12ap and 

dQ = Q x u x + Quu u + Q P oj p + A(0) 

= Qx u x + Qu co u + Q P uj p + v X x + Qiyu ~^x) + P{^xu - wx) 
- WQnu + P 2 {v uv - 2nxu) - P'^uu, 
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/ du 


— pdx\ 






1 dp 


— qdx \ 





a 4 / 




dx j 



we see that it is possible to set 

X, U, P, Q -> (3.3) 

as their exterior derivatives involve the linearly independent Maurer-Cartan forms (i, 
v, fx, vxx- The result is the system of equations 

= u x + H, = cj u + v, = w p + z/ x , = Q P uj p + Q v uj u + Q x uj x + v X x, 

which can be solved for the partially normalized Maurer-Cartan forms: 

H = -uj x , u = -cj u , u x = -uj p , vxx = -(Qpoj p + Quoj u + Q x u x ). 

Remark: after substituting (I3.3D in (I2.5D and (I2.8P and normalizing the pseudo-group 
jets, we recover the usual G-structure formulation of the equivalence problem |15j . 



(3.4) 



where 

ai = U u , a 2 = D J* y" , a 3 = X u , a 4 = D X (X) = pX u + X x . 

Continuing the normalization procedure, order by order, let H C J'^- 00 ^ be the sub- 
manifold defined by 

x = u = p = 0, (3-5) 
qij0 = qiji = 0, i,j>0, 

qoj2 = qij2 = <?oj3 = qijs = 0, j > o, 

and let 

£ = t~ 1 (E)c£ {co) 
denote the lift of H; that is, £ is defined by the equations 

X = U = P = 0, 

QijO = Qiji = 0, i,j>0, (3.6) 
Qoj2 = Qlj2 = Qoj3 = Qlj3 = 0, j > o. 

Let H C SL3IR be the 5-dimensional subgroup 







«2 








!(• 




! ) 


: ai6 2 C3 






C2 


c 3/ 





of fractional linear transformations that preserve the origin (x,u) = (0,0). The proof 
of the following Propositions is presented in Appendix [Bj 

Proposition 3.1. The submanifold H C J'( 00 ^ is a global cross-section for the action 
of Diff(M 2 ) on J(P°). 

Proposition 3.2. We have H = Q^l that is, H is the subgroup that preserves 3. 
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In light of Proposition 13. 1\ we can utilize S as a partially normalizing cross-section 
for the point-equivalence problem. Given a smoothly defined ODE q = f(x,u,p), 
let df. £ f -> N 3 denote the pullback of a : £ -> jt 00 ) via /(°°) : iV 3 -> J^ 00 ). By 

Proposition 13.21 this pullback is a reduction of structure from a f : cf° — > J^ 00 ) to the 

principal i/-bundle, 07 : £/ — >■ N 3 . 

In the sequel, we use the tilde decoration to denote the pullback to £f, and refer to 
the quantities 

£ f 

as universal invariants. The universal invariants are, in fact, the components of ff : £f —> 
E; the latter obtained by imposing the normalizations (|3.6p . As such, the universal in- 
variants are functions of 3 + 5 = 8 variables and are iJ-equivariant with respect to the 
restricted left action (I2.ip , 

As is shown in Appendix [B] the following one-forms 

fj,, v, i/ x , "xx, fj,ij,i + j>2, Vij,i + j>3, 



are normalized by (|3.6p . In particular, applying (|3.6p to the universal recurrence for- 
mulas for 

dQpk, dQpkx, dQpkjj, fe = 0,l,2,3, 
yields following relations: 

fi = —£j x , v = —Cj u , vx = —& p , 

"x 2 = "X 3 = ^x 2 u = fix 4 = 0, fix 2 = Zvxu, 

1. „ 1~ „ „ ~ (3-7) 

fixu = fljji = -Q P 40J P , 3fix 2 u = 6v X u 2 = Qp 2 x 2 & x , 

2 b 

6fixu 2 = = Qp3x 2 & x + Qxp 4 & P i 6fiu 3 = —Qp 4 &xu + Qp 4 u ^ v ■ 

Substituting the above relations into (|2.13p and (j3.2[) yields the structure equations for 
£f, namely: 

doj x = fix A u) x + fiu A u> u , 

du) u = %A(i 11 + w 1 A oj p , 

doj p = vjjx A Cj u + (% - fix) A cD p , 

d/ix = -2 %x Aii 1 - -%[/ A w 1 * — fiu A cD p , 

= —zf>uu A w z + (/2x - A /it/ + ^Q P 4W M A U) p , ( 3 - 8 ) 
2 b 

d% = —Dux A w 1 + flu A lD p — %{/ A a)", 

d%(7 = 2%x A fiu + vuu A % + \^Qpix & u A w p + ^Qpaxa w" A 

1 1 ~ 

d^ax = z^c/x f\fix- -vuu Aw p + -Qp2x 2 ^ u A w 1 . 



We now define the reduced rank sequence 

/ 



£•„ := rankf ( , n) , n > 0. (3.9) 
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Proposition 3.3. The following are equivalent: 

(i) Q — f( x i u iP) is a regular ODE; 

(ii) for every n, the reduced rank g n is constant; 

(hi) the invariant coframe oj x , & u , oj p , £t x , fiu, %> ^xu> ^uu on £f is fully regular [T5l 
Dehnition 8.14]. 

Proof. The equivalence of (i) and (ii) follows from the equivariance property (|2.2|) of the 
target projection. The equivalence of (ii) and (iii) follows from the structure equations 
(13.81 ) and from the universal recurrence relations ()2.12bj) . □ 

In light of the above remarks and Proposition 13.31 we observe that 

Qo = Qi = 02 = Q3 = 0, and g 4 < g 5 < g e < . . . < 8. 

It follows that the reduced rank sequence stabilizes at a sufficiently high order n. Indeed, 
the following is true. 

Proposition 3.4. The IC order can be characterized as the smallest integer n > 4 
such that g n -\ = g n . 

Proof. See Proposition 8.18 of [15]. □ 



3.3 The fundamental branching 



The universal reduction (I3.6P leads to the normalization of all the Maurer-Cartan forms 
(|2.6p except for 

fix, fiu, %, thru, vxu- (3.10) 
To proceed further, the value of the universal invariants 

Qpk+4, V j X i, Q P 3 U j X i+ 2 , Qp 2 u^x i + 2 i i,j,k>0, (3-11) 

(recall that these are restrictions of the lifted invariants Qijk to the universal normalizing 
cross-section (|3.6j) ) must be analyzed in more details. Up to order 6, the non-trivial 
universal invariants (|3.1ip are 



n = 4 
n = 5 
n = 6 



Qp4, Qp2 X 2, 

Q P 5, Qpijj, Qpixi Qp^x 2 , Qp 2 ux 2 , Qp 2 x 3 , , , 

Qpe, Qpsjj, Qpsx, Qp^u 2 , Qp^uxi Qp^x 2 , Qp s ux 2 i 



Q P 3 X 3, Qp2u2 X 2, Qp2 UX 3, Q p2 X A . 

Writing the recurrence relations for the invariants (|3.12j) of order < 5 we obtain 
dQpi = Q P s Cj p + Q P 4ijCj u + Q P 4 X Cf + Q P i{2fi x - 3%), 

dQp2 X 2 = Qp3 X 2 0J P + Qp2jJ X 2 Ld U + Qp2 X i U) X — Qp2 X 2 {y\j + 2fix), 

dQ P 5 = Q P e Cj p + Qp5ijCj u + Q P 5 X Cj x + hQ P i flu + Qps(3/ix - 4%), 

dQp4 X = {Qp5 X + QpAjj) & P + Qp^UX & U + Qp4 X 2 UJ X + Qp4 VJJX 

+ Qp+xifix - 3%), 
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dQpiu = Qpsjj £j p + Qp4jj2 u) u + Qpijjx & x ~ "2Qpi — Q P s vjjx ~ Qp 4 x fiu 

+ Q P 4 U {2fi x -4%), (3.13) 
dQp-Ax 2 = Qp 4 x 2 & p + Qp 3 ux 2 & + {Qp 3 x s ~ 2Qp 2 ux 2 ) & x — Qp2x 2 fiu 

- Qpsx 2 {^u + fix), 

dQ P 2jjx 2 = Qp3jjx 2 & v + Q P 2jj2x 2 ii" + Qp2 UX i & x — 2Qp 2 x 2 ^uu — Qp^x 2 &ux 

- Qp 2 X3 fiu - 1Qp 2 UX 2 (vu + fix), 

dQp2 X 3 = (Qp3 X 3 — Qp2 UX 2)Cj P + Qp2 UX 3 u" + Q p2 X 4 — 5Qp2 X 2 i>ux 

- Qp 2 x s ( i> u + 3fix)- 

At this juncture, the equivalence problem splits into 4 branches depending on the values 
of the fourth-order universal invariants: 

I) Q P 4 = and Q P 2 X 2 = 0, III) Q P 4 = and Q P 2 X 2 ^ 0, 

II) Q P i ^ and Q P 2 X 2 = 0, IV) Q pi ^ and Q P 2 X 2 ^ 0. 

Branch I corresponds to the equivalence class of linearizable differential equations 
discussed in the introduction. For this class of equations we have Q4 = 0, and so the 
IC order equals 4. 

As for case IV, we see from the recurrence relations (|3. 13[) that the Maurer-Cartan 
forms (|3.10p can be normalized by setting 

Q P 4, Q P 2 X 2 — > 1, Q P 5,Qp4 U ,Qp4 X — > 0. 

This uses up all of the remaining //-freedom and produces a genuine moving frame. The 
algebra of absolute differential invariants is then generated by the remaining invariants 
(|3.12p of order 5 and 6. The reduced ranks are 

£?4 = 2, £? 5 > 5. 

Hence, the "worst-case scenario", as far as the IC order is concerned, is 

{h, Qs, h, h, h) = (5, 6, 7, 8, 8); 

and as a consequence, the highest IC order achievable is 9. This bound will be attained 
if, post-normalization, the remaining 5th order invariants are constant, and the signa- 
ture manifold is parametrized by three invariants of order 6, 7, 8, respectively, with the 
higher order invariants obtained by differentiating the order 6 invariant. We do not 
push the analysis further as we will show that cases II and III contain equations with 
invariant classification order equal to 10. Since our goal is to find the branch(es) with 
highest classification order, we will focus on those branches. Indeed, in the sequel we 
consider case III in detail, as branches II and III are dual to each other. This duality 
was first observed by Cartan in his study of projective connections, [I]. For complete- 
ness, the duality among second-order ordinary differential equations is presented in 
Appendix [Al 

3.4 Case III 

From now on, we assume that Qp2x2 ^ and Q P 4 = 0. With these assumptions, we 
will show that the reduced rank sequence obeys 

Qa = 1, 05 = 4, £?6 > 5. 
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Hence, for this class of ODEs there exists a genuine moving frame formulated in term 
of 6th order jets. For class III equations, the "worst case scenario" is the rank sequence 

(f?4, 65, Q6, 67,08, 89, Qw) = (1,4,5,6,7,8,8); 

which makes an IC order of 10 a possibility. 

By Proposition 13.51 below, the class III equivalence problem has two branches, 
which we label as III.l and III. 2. For sub-case III.l we will show that qq = 5 implies 
Q7 = 5, which means that the IC order is 7. The other possibility is that qq > 6, but 
this means that the IC order is < 9. Hence, sub-case III.l can be ruled out. 

Finally, for sub-case III. 2 we will show that there is essentially one type of configu- 
ration of invariant values that gives a rank sequence of (1,4, 5, 6, 7, 8, 8). We will derive 
this configuration, and in the subsequent section integrate the corresponding structure 
equations. 

Under the non-degeneracy assumption Q P 2 X 2 ^ it is possible to normalize 

Q P 2 X 2 -> 1, Qpi X 2 ,Qp 2 ux 2 ,Qp 2 x^ 0, (3.14) 

which consequently normalizes the Maurer-Cartan forms Djj, v\j, i>uu, ^ux to cer- 
tain linear combinations of oj x ,ui u ,u} p ,flx- Henceforth, to avoid confusion, we use the 
"check" Qijk decoration to indicate the invariants and one-forms obtained via additional 
normalization of the universal invariants as per (|3. 14[) . Furthermore, the recurrence re- 
lation for Qp4 = 0, forces the following fifth-order invariants: 

Q P 5 = Qpijj = Qp±x = (3.15) 

to be identically equal to zero. Combining f|3. 14j) and (|3. 15[) we conclude all universal 
5th order invariants can be normalized to a constant, and that to normalize the re- 
maining Maurer-Cartan form fix we must consider universal invariants of order 6. At 
order 6, the constraints (I3.15P force the invariants 

Qpe = Qpsjj = Qp$x = Qp^u 2 = Qp 4 ux = Qp^x 2 = (3.16) 

to be identically zero, and more generally, 

Q P A +iU j X k = 0, i,j, k > 0. (3.17) 

Hence from (|3. 16[) we conclude that the remaining non-constant universal invariants of 
order 6 are 

Qp^ux 2 , Qp^xz, Qp 2 u 2 x 2 , Qp 2 uxz, Qp 2 x 4 - 

Proposition 3.5. The invariants Qp2 X i and Qps X i cannot simultaneously be equal 
to zero. 

Proof. Considering the recurrence relations for Qp3 X 3, Qp 2 ux 3 , Qp 2 x 4 we have 

dQ P 3 X 3 = -^Qpiux 2 & p + {Qp^ux 3 ~ ■^Qp 2 u 2 x 2 ^Ju u + (^Qpi X i - -^Qp 2 ux'^ju x 
+ 2Qp3x3 fix, 

dQp 2 ux 2 = (Qpzuxz - Qp 2 u 2 x 2 - -^Qpsx^j^ P + {Qp 2 u 2 xz - Qpzux 2 Qp 2 x* 
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- -^Qp»x^ Qp 2 ux^j^ u + (Qp 2 ux± - g - gQpaxs Qp2x4^w ;I + Q P 2 UX i Ax, 

dQp-2xi = <5p3jf4 w p + Qp2jjx 4 & u + Qp 2 x 5 & ~ 2Qp?x 4 fix- 

Assuming Q P 2 X 4 = Qpsx 3 = 0, the syzygy 

9 „ 

D x {Qp3 X s) + ^Qpaf/xs = Qpix 4 = Dp(Q P 2 X i) 

forces Qp2ij X 3 = 0, which when combined with the syzygy 

5 6 - 

D x (Q P 2 UX a) + - + -<5p3x3Qp 2 x 4 = Qp 2 ux* = Du(Q P 2 X 4), 
o 5 

leads to the contradiction = 5/6. □ 
By virtue of Proposition 13.51 two sub-cases must be considered: 

III.l) Qp3 X 3 # 0, III.2) Qp2 X 4 0. 

3.4.1 Sub-case III.l 

Assuming Qpz x z ^ 0, we can set 

Qp3 X 3 —7" 1 

and normalize the Maurer-Cartan form /ix- After normalization, the remaining sixth- 
order invariants 

Qpzux 2 , Qp 2 u 2 x 2 -> Qp 2 ux s > Qp 2 x 4 , (3.18) 

are genuine invariants in the sense that they do not depend on pseudo- group parameters. 
In an attempt to minimize the rank, we assume that the functions (|3.18p are constant: 

Qp s ux 2 = Cij Qp 2 u 2 x 2 = C2, Qp 2 ux 3 = C3, Qp 2 x± = Ca. (3.19) 

Combining ()3.16j) with (|3. 19p . it follows from a careful analysis of the recurrence rela- 
tions that all seventh-order invariants are constant which in turns forces all higher-order 
invariants to also be constant. On the other hand, if the invariants (|3. 18[) are not con- 
stant, the IC order is < 9. 

3.4.2 Sub-case III.2 

We now assume that Qp2 X 4 ^ 0, and set 

Q P 2 X 4 -> 1 (3.20) 

to normalize fi x and obtain a genuine moving frame. From now on, for the sake 
of notational convenience, we omit writing the check decoration and simply use Qijk 
to denote the absolute differential invariants obtained by normalizing the universal 
invariants using f)3. 14j) and (|3.20p . Likewise, the invariant coframe on N 3 will be written 
simply as uj x , uj u , uj p and the dual derivative operators as D x ,Djj,Dp. 
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At order 6, we are left with the absolute differential invariants 

QpZUX 2 , Qp3 X 3, Qp2 V 2 X 2, Qp2 UX 3- (3-21) 

Once more, in an attempt to minimize the rank, we assume that the invariants (|3.2ip 
are constant: 

Qp^ux 2 = C\, Qpsx'> = C2, Qp2 U 2 X 'i = C3, Q P 2 UX 3 = C 4 . (3.22) 

The seventh-order invariant Q P 2 X s plays an important role in the following considera- 
tions. To single out this invariant, let us set 

I7 = Qp2 X 5. 

Considering the recurrence relations of the invariants (|3,22p . we obtain a collection of 
constraints on the seventh-order invariants: 

= d Cl = 5 -^f^ ^ + (q p3u2x2 _ ClC2 + SCiQpw ) u u 

+ (QpWx 3 2C3 H Ij-^ju'- 1 , 

= dC 2 = (c 2 Q P z X 4 + ^)^ p + (Qpzux* + C 2 Q P 2 UXi - ^ V 

+ [Qp^ +I7C2 - (3.23) 

Qp-i V 2 X 2 + 2C 3 Qp3 X 4 \bJ P + ( Qp2 U 3 X 2 + 2C 3 (5p2 l/X 4 

- V« + (q p2u2x3 + 2/ 7 C 3 - 2C 2 C 4 - 2 -^y 
= dC 4 = {qp^ux^ - C3 H — 4< ^f X s~) u,P + ^Q^ 2 ^ 2 ^ 3 ^ x 

"—J" +^^ + -2 r~6j u • 

If the invariant Ij is constant, (|3.17p and (|3.23p imply that all seventh-order invariants 
are constant. Similarly, all higher order invariants are constant. On the other hand, 
when I j is a non-constant invariant, the constraints (|3.23p yield 

C\ = C 2 = C3 = C4 = 0. 

Which in turn, implies, together with (|3.17p . that all seventh-order invariants are iden- 
tically equal to zero except for 

5 

Qp2 X s = I7 and Q P 2 UX a = -. 

6 

Taking the exterior derivative of Ij we obtain 

dI 7 = Izu*-- A I 7 u u + \u*>, (3.24) 

4 D 



CO 
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where 1$ = Dx{Ij) is the only new (functionally independent) invariant of order 8. 
Then, differentiating 1$ with respect to Dx we find the only new invariant of order 9: 

h = D x {h) = D\{I 7 ). 

Generically, the invariants I 7 ,I$,Ig, are functionally independent, and the structure of 
the invariant signature manifold is completely determined by the functional relation 

I w = D x {h) = 4>(I 7 ,h,h). 

Modulo duality, all branches of the equivalence problem have now been considered, and 
we can safely conclude that 10 is an upper bound on the IC order. 

4 The maximal IC order class 
4.1 Abstract existence 

To terminate the proof of Theorem 11,31 we must show that there exists a class of 
differential equations satisfying the invariant constraints imposed in sub-case III. 2. To 
do so, we need the structure equations of the invariant one-forms u' x , ui u , d p . 

These equations are obtained symbolically by substituting the Maurer-Cartan form 
normalizations 

fi = —co x , v = —u u , fx = — w p , 

5 1 

-2fx x = v v = -I 7 ui x - - u u , nu = vxx = 0, u ux = -u x , 

o 5 

into the structure equations (|3.2p . The result is 

5 

dw x = — U) u A uj x , 
12 

duj u = ui x A up + Ij uj u A ui x , (4.1) 

dio p = -I 7 oj p A co x + - co p A uj u + - co x A to 11 . 
2 4 5 

If the invariant I7 were constant, we could apply Cartan's Integration Theorem, [3], 
to conclude the existence of differential equations solving the integration problem. But 
since this is not the case we must use the following generalization of Cartan's integration 
theorem [2]. 

Theorem 4.1. Let ui 1 , . . . , ui p be a coframe with structure equations 

i<j<k<p 

such that the structure coefficients are functions of I a , a = 1, . . . , s and 
dP = ( F ?(I b ) + £ AUl b ) J a ) u\ a, b = 1, . . . , s. 

i=l ^ a=l ' 

Assuming 
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• the identity d 2 = holds, 

• the functions C* fc , F®, Af a are real analytic, 

• the tableau A(I b ) = (A? a (I b )) has rank r and is involutive with Cartan characters 
si > S2 > • • • > s q > s q+ i = for all values of (I b ), 

modulo a diffeomorphism, the general real-solution exists and depends on s q functions 
of q variables. Moreover, (I a ) and (J a ) can be arbitrarily specified at a point. 

Applying Theorem 14.11 to (|3.24p and (j4. 1 j) we conclude that there exists a family of 
second-order ordinary differential equations depending on one function of one variable 
that solves the integration problem. Hence, by duality we conclude the there are two 
families of equations, both depending on one function of one variable, that achieve the 
maximal classification order of Theorem 11.31 



4.2 Explicit integration 

In this section, we integrate the structure equations (14. lj) to obtain an explicit repre- 
sentation of one of the two families the differential equations satisfying the maximal 
classification order of Theorem 11.31 

Proposition 4.2. Let (a, /3, 7) be a local coordinate system on M 3 \ {a = 0}, then the 
one-forms 

r , ,/ 12 da n , „ d/3 12 1 7 da n . , , t 

w x = ad7, w M = - /3c?7, ^ p = — + — -— + r(a, /3, 7) dj, (4.2) 

5 a a 5 a 

where r(a, /S, 7) is a solution of the of linear partial differential equations 

3 12I 7)7 18/ 7 2 3/ 7 /3 12 3/ 7 5/3 

+ = — 1 = — , lp = — — , {4.6) 

a ba 5 a 25 2 4a 

satisfy the structure equations (|4.ip . 

Substituting (14.21) into (I3.24p we obtain the differential equations 

j 5 j 7i7 5/3/ 7 sr 

-'7, / 9 — 7—, ai 7a — — i 7 , is — -. — • 

ba a 4a ba 

Integrating the first two equations, we deduce that 

I 7 = ¥ + k M ,4.4) 
ba a 

where ^(7) is an arbitrary (analytic) function, while the third equation defines Is in 
terms of I 7 and T. The coframe (14. 2D is not uniquely defined. The degree of freedom is 
given by the infinite-dimensional Lie pseudo-group 

- r ^ - a -r P l2a " 
7 = cr(7, a = — , p = — - — 

a' a' b{a')' L 

-_ T 1217a" fia" 12a 1 " 2\{a"f ^ 
r " a 7 + 5(cj') 2 + ^(V) 1 + 5a(a') 2 ~ Sa^') 3 ' 
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where 7 = 17(7) is a local diffeomorphism of the real line. Under the pseudo-group 
action (|4.5p . the invariant (|4.4p transforms according to 



5/3 2a" a'h 
6a aa' a ' 



I7 = H , with h = h°o~. 



By choosing (7(7) such that 



2a" = (a') 2 ^, 



we can assume ^1(7) = in (|4.4p . Doing so and solving the differential equations (|4.3p 
we find that 

r (a , ft7) = -^ + £W, 

where 5(7) is an arbitrary (analytic) function which cannot be removed by some change 
of variables. 

Proposition 4.3. The one-forms 

u* = a dl , ^ = ^-^7, ^ = ^ + 2 4da + ( 9 M_^) dl (4.6) 
5 a a ar \ a° 25 / 

satisfy the structure equations (|4.ip with 

I(a,P,i) = ^-. (4.7) 



Proposition 4.4. Let uj x , co u , lo p be given by (14. 6p . Then there exist functions 01, 02, 
03, 04 and a change of variables 

a = a(x,u,p), ft = P(x,u,p), j = y(x,u,p) (4.8) 

such that (13.41) holds. 



Proof. The functions ai, 02, 03, 04 and the change of variables (|4.8p are not unique. 
Assuming u > 0, we can choose 

6 

ai = — ai = CI3 = 0, 04 = a, 
oa* 



and 



Then 



a = y/u, p = — , 7 = x. 
5u 



oj x = u l l 2 dx, uj u = — [du - pdx], uj p 



5u 5u 3 / 2 



dp-(^-^)dx 
F V 10 6 



• (4-9) 
□ 



We deduce from co p in (14. 9h that the second-order ordinary differential equations 

u 2 5g(x) 



10 6 



(4.10) 



is one of the 2 families of differential equations with maximal classification order, pro- 
vided I7,l8,lg are functionally independent. A straightforward scaling transformation 
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gives the form shown in (|1.7|) , We note that the coordinate expressions for the dual 
equations to (|4.10p are very difficult to obtain. 

Finally, for the differential equations (|4.10|) . the coordinate expressions of the in- 
variant (|4.7p and the derivatives dual to oj x ,lj u ,u p are 

(4.11) 



' u 3 / 2 
and 



D x 



ii 



1/2 



d d ( u 2 5q(x) \ d 



dx du V 10 6 J dp 



hu d 5ti 3 / 2 d 

Du = ~pr ~r\ j Dp 



6 du 1 6 9p' 



respectively. Hence, differentiating (|4.1ip twice with respect to Dx we deduce that the 
invariant signature manifold can be parametrized by the invariants 

/7 = ^72' h = ^-, h= 5f , = 256 [D x {g(x) 1 )) . 

The above invariants are independent if and only if Ig is non-constant, which gives 
the inequality in (jl.7p . If this inequality holds, then the invariant classification of an 
equation in this class is completely determined by the functional relationship between 
Ig and the invariant 

g{x)g"(x) 

A Cartan's duality for second-order ODEs 

In his study of projective connections, [1], Cartan mentions that there exists a notion 
of duality among second-order ordinary differential equations. Modern accounts can be 
found in [U[T3], but for completeness, we summarize the construction below. 

For second-order ordinary differential equations, a solution to an initial- value prob- 
lem may be represented by a two-parameter family 

$(x,u,x,u) = 0, (A.l) 

where the parameters x, u correspond to the initial values 

x = it(0), u = u x (0). 

Differentiating (jA.ip twice with respect to x we obtain the equations 

$ = Sj. + u x $ u = § xx + 2u x $ xu + ul^uu + u xx <& u = 0, (A.2) 

which leads to the second-order differential equation 

u xx = q(x,u,u x ) (A.3) 

when the parameters x and u are eliminated from (|A.2p . The dual equation to (|A.3D 
is obtained from (jA.ip by interchanging the roles of (x,u) and (x, u). In other words, 
let (x, u) be parameters and u = u(x) a function of the independent variable x. Then, 
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differentiating (|A.ip with respect to x twice and eliminating (x, u) from the equations 
obtained yields the dual equation 

u-xx = q{x,u,Ux). (A.4) 

Let us now consider the consequences of the contact transformation that sends (|A.3P 
to (|A.4p on the point equivalence problem. For the equations 

&(x, u, x, u) = 0, & x (x,u,x,u) + p& u (x,u,x, u) = 0, 

to determine x, u, p in terms of x, u, p we must impose 

/ A = det 3^ ^ $ x . s = -Sv^v^s + + P^u + PP$t«z), 

which in particular requires 

$ u / 0, $^ 0. 

Let 

9 = du — palx, 9 = du — p dx, 9\ = dp — q dx, 9\ = dp — q dx. 
Then, taking the exterior derivative of (|A.1|) we deduce that 

9 = -^9. 

On the other hand, the exterior derivative of & x +p& u = yields 

A 



dx mod 9, 



while the exterior derivative of $x + P &u — gives 



dx = -Jt-Ji 0i mod 0. 
A 



In matrix form, 



-ai 



03 04 



dx I = I a2 01/04 I I 0i I , where oi = — -, 04 



We note that, up to a sign in the first entry, the 3x3 matrix made of the functions 
Q>1, 0,2, 03, 04 is an element of the structure group (|3.4p . Hence, under the contact 
transformation (x,u,p) (x, u, p), the lifted coframe f|2 . 5f) and (j2.8H undergoes the 
transformation 

lo x o w 17 , w p <-> w*, w" o -u*. (A.5) 

To understand the duality between cases II and III of the equivalence problem in 
Section 13.31 we consider the structure equations (|3.8|) obtained once the normalizations 
(|3.6p are done. Focusing on the structure equations 

1 ~ 1 ~ 

dfiu = 1- -Qp4 li™ A Co p , dvux = h -Qp2 X 2 w u A 

6 6 

we observe that under the coframe transformation (IA.5h the role of the universal in- 
variants Qp4, Qp2x2 is interchanged, which is exactly what happens when switching 
between cases II and III. 



21 



B Proofs of Propositions 13.11 and 13.2 



We start by proving Proposition 13. 11 Let us first consider the local transversality. First, 
we observe that dX, dll, dP, dQijk, fiij, v^a and oj x ,lo u , uj p , 'dijkilHjt v ij are two choices of 
right-invariant coframes on £(°°\ Then, we note that the left-action of Diff(IR 2 ) in (|2.1|) 
generates the source distribution of a : — >■ i y( 00 ) given by 

keic{dx, du, dp, Oijk} = kerfo/ 1 , w", u p , i?^}, 

while the right-action generates the target distribution of r: — >• J'^°°\ given by 
kei{dX, dU, dP, dQijk}- Hence, our goal is to show that the tangent space to £ , defined 
by equations f|3.6[) . is transverse to the source distribution on £^°°\ 
According to (12. 9p . 

<f>(x, u,p, q) = r] xx + q(rj u - 26) + p(2r] xu - £ xx ) - 3pq£ u + p 2 (r] uu - 2£ xu ) - p 3 £ uu - 
Writing 

& ] 0<i + 3 <n}, 

Vij = Vx^ui , V (n) = {Vij : < i + j < n}, 

the prolongation formula ()2.10j) entails 



< 


jOOO 


= Vxx 




mod p, C^,?/ 1 '; 




b ooi 


= 2t\ xu 




mod p, j 


<i 


,,002 


= 2t] uu 




mod p, ^'jV , 


S^j-2,2 _ ^IJ 


-3,3 


= 6??oi 




mod p,^ -1 ^ -1 ), j > 3; 




-2,2 


= Qmj 




mod p, , Tyi' , j > 2; 


<t 


2J0 


= Wj 




inodp,C (i+i ~ 1) ,?7 (<+i ~ 1) ) » > 2,j > 0; 




-2,3 


= -6£oj 




mod p,£ (j_1) ,7? (i_1) , J > 2 ; 




-2,3 


EE -6^y 




mod p, £^ , rft' , j > 2; 




-1,2 


= 36, 




mod p,£ (i+1) ,r/ (i+1) , j > 1; 


p-2j,l _ 2 0i-3j 


fl,0 






modp,C (i+J '" 1) ,?7 (i+i ~ 1) ) » > 3,j > 2. 



By the universal recurrence formulas ()2.12p . modulo P, co x , uj u , co p , $ijk we have 



dX 


= fj,; 








dU 


= v\ 








dP 


= vx\ 








dQooo 


= v X x 




mod yP^ 




dQooi 


= 2u X u 


— l^xx 


mod /i^ 


u m. 


dQoo2 


= 2vuu 


- tyxu 


mod //^ 




3dQo,j-2,2 — %dQi t j-3 t 3 


= 6v j 




mod fjy~ 


j>3 


4dQoj-i,i — dQ\,j-2,2 


= 6 ^ij 




mod fj,W 


i/W, j>2; 
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dQi-2,jO = v ij 
dQo,j-2,3 = — 6Moj 
dQl,j-2,3 = — 

dQo,j,i — dQij-ip = 3/Lt2j 
dQi-2,j,i — 2<i(5i-3,j+i,o = /%' 

It follows that uj x ,uj u ,uj p ,vjj, [ix-, i v xu i l/ ujj-,'&ijk form a basis of T*£. As a conse- 
quence, since, u u , lo p , are linearly independent on £, the latter is transverse to 
the source distribution. 

Now we prove the surjectivity. Let J = {(i, j, k) : k < 2 or (i < 2 and k < 4)} be 
the indicated set of indices. We have to show the consistency of the following equations: 

P(p,X x ,X u ,U x ,U u ) = 

which is equivalent to 

pU u + U X = 0; 

and 

Since Q is transitive on N 3 = J" 1 (R, R), no generality is lost if we set p = above; that 
is, it suffices to demonstrate the consistency of the equations 

= Q t ^q { ^ +k \x^ +k+2 \ut J+k+2) ), (i,j,k) e J. 

Furthermore, the above equations can be greatly simplified by restricting to the sub- 
pseudo-group of Q defined by 

X X = U U = 1, U X = X U = 0. (B.2) 

At p = the transformation law for Qijk, k < 4 is, modulo (|B.2|) and X^ i+ ^ +1 \ L7( i +J+ 1 ), 
affine to leading order 

QijO = QijO + Ui + 2,j, 
Qijl = + 2C/j + i J+ i — Xj + 2J, 

Qij2 = Qij2 + 2Uij + 2 — 4Xj + ij + i, 
Qij'3 = 9ij3 ~~ 6Xjj + 2. 

We now impose relations for which (i, j, k) £ J. Setting Qififi = Qi,j+i,o = Q%+i,j,i = 
fixes U i+ 2fi,U i+ 2,j + i, X i+3J . Setting Q 0j3 = Q lj3 = fixes X j +2 , Xij +2 . Fixing 
Qoj+i,2 = Qij+i,2 = determines t/oj+3, Lqj+3, and then setting Qo,j+2,i = specifies 
-^2,?+2 ■ Consequently, the remaining relations are 



Qooi 


= q oi + 2f7n 


- X20, 


Qon 


= gon + 2f/i2 


— X21, 


Ql02 


= q W2 + 2U 12 


-4X21 


Q002 


= Q002 + 2?7o2 


-4X n 


Q012 


= Q012 + 2Uq 3 


-4X 12 



mod ^i- 1 ),!/^'- 1 ), » > 2,j > 0; 
mod j > 2; 

mod y^\u^\ j > 2; 
mod A t&' +1 ),i/(j +1 ), j > 1; 
mod ^ i+ i +1 \u( i+j+1 \ i > 3, j > 2. 
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Setting the left-hand side to zero, these relations fix X 12 ,X 21 , ?7i2,-^20)-^ii leaving U~n 
and U02 as free variables. 

We now consider the proof of Proposition 13.21 Let 



£e = n r~ 1 (H) C £ 

be the subgroupoid of transformations that preserve 3, and let £h be the subgroupoid 
corresponding to the actions of H on J'(°°\ It is straightforward to check that H pre- 
serves 3; that is, C £=• Equations (IB.lj) also demonstrate that vy , /xjf , fj,u , ^xu, v vv 
form a one- form basis for the source fibres of a : £~, — > 3. Since is 5-dimensional, it 
follows that Eh = £z by dimensional exhaustion, as was to be shown. 
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